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Abstract. The dynamic patterning of the plant hormone auxin and its efflux facilitator the PIN protein 
are the key regulator for the spatial and temporal organization of plant development. In particular auxin 
induces the polar localization of its own efflux facilitator. Due to this positive feedback auxin flow is 
directed and patterns of auxin and PIN arise. During the earliest stage of vein initiation in leaves auxin 
accumulates in a single cell in a rim of epidermal cells from which it flows into the ground meristem tissue 
of the leaf blade. There the localized auxin supply yields the successive polarization of PIN distribution 
along a strand of cells. We model the auxin and PIN dynamics within cells with a minimal canalization 
model. Solving the model analytically we uncover an excitable polarization front that triggers a polar 
distribution of PIN proteins in cells. As polarization fronts may extend to opposing directions from their 
initiation site we suggest a possible resolution to the puzzling occurrence of bipolar cells, such we offer 
an explanation for the development of closed, looped veins. Employing non-linear analysis we identify the 
role of the contributing microscopic processes during polarization. Furthermore, we deduce quantitative 
predictions on polarization fronts establishing a route to determine the up to now largely unknown kinetic 
rates of auxin and PIN dynamics. 

PACS. 87.17.Pq Morphogenesis - 87.10. Ca Analytical theories - 87.18.Hf Spatiotemporal pattern forma- 
tion in cellular populations - 82.40.Bj Oscillations, chaos, and bifurcations 
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1 Introduction 

The polar transport of the plant hormone auxin is the key 
regulator of many processes in the spatial and temporal 
organization of development and growth of plants. As the 
indole-3-acetic acid, in short auxin, induces the polar lo- 
calization of its own efflux facilitator, a member of the 
family of PIN proteins, a variety of auxin and PIN pat- 
terns arise [T]. Those distributions change dynamically as 
plants orient in response to environmental stimuli denoted 
tropism 2,3. During the morphogenesis of plants PIN and 
auxin rearrangements lie at the heart of organ positioning 
via phyllotaxis @] and vein patterning in leaves [5J. 

The notion that auxin is transported in a polar, di- 
rected, manner inspired researchers since its discovery by 
Went in 1933 [6 . Early works already suggested the partic- 
ipation of a polar localized efflux carrier in the transport of 
auxin 7,8,9, well before its discovery in the form of mem- 
brane bound PIN proteins a decade ago [TU]. Since then 
numerous experiments confirmed that PIN proteins facili- 
tate the efflux of auxin from cells in plants [TTl[2l [T2l[T3T ll] , 
yeast and mammalian cells, which had been supplied with 
auxin and PIN [TJ]. A feedback between auxin and its ef- 
flux facilitator localization was proposed by Sachs in his 
canalization hypothesis p~5j, later formalized by Mitchi- 



son [16, 17j. Canalization predicts a feedback of auxin flow 
between neighboring cells on the amount of efflux facili- 
tators favoring the direction of auxin flow. Experiments 
confirmed a definite feedback between auxin and PIN dis- 
tribution [P8"lll9| . the cause of which is reported to lie in 
auxin affecting the clathrin dependent endocytotic cycling 
of PIN [20 1 121 ] . Late investigations also identified biochem- 
ical processes taking part in the PIN localization in re- 
sponse to auxin, see Ref. |22] for a review. 

A variety of microscopic models for the dynamics of 
auxin and PIN proteins have been developed to describe 
their patterns during phyllotaxis [33] and leaf vein for- 
mation [24] . see [15] for a review. Extensive simulations 
of these microscopic models describe qualitative aspects 
of plant development. However, the role of the underlying 
biological processes and their kinetic rates still remain elu- 
sive to a large extent. Quantitative predictions based on 
analytical solutions of the microscopic equations in a sim- 
ple scenario might on the one hand help to estimate kinetic 
parameters and on the other hand give insight into the im- 
pact of certain processes. A scenario amenable to such an 
investigation is the polarization of in this particular case 
PIN1 distribution due to auxin flow in the earliest stage 
of vein formation [55] , see fig. [I] Vein initiation itself takes 
place in the ground meristem tissue of leaf primordia. The 
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Fig. 1. Illustration of the dynamics of auxin and its efflux facilitator PIN. (A) Schematic drawing of vein initiation in a leaf 
primordium. Auxin (blue) accumulates at convergence points in the outer epidermal layer (rim), from which it is transported 
into the ground meristem. Due to this inflow of auxin ground meristem cells become polarized in their PIN distribution (orange). 
(B) Polarization of PIN distribution in a strand of cell due to auxin inflow from the left, indicating details of auxin and PIN 
dynamics. The weak acid auxin accumulates in the interior of a plant cell due to a gradient in pH. In the inner cell the charged 
anion is trapped and can only be transported outwards by help of efflux facilitators in the form of PIN proteins. The auxin 
transport from cell to cell has an efficiency ba- Auxin synthesis sa and degradation cLa takes place in the inner cell. PIN proteins 
cycle between the bulk and the cell membrane by basal attachment sp and detachment rates dp. In addition, positive auxin net 
flow is modeled to feed back on these rates increasing the PIN attachment by gp. Along the strand of cells the color shading 
indicates relative concentration of auxin and PIN on either membrane. 



positions of vein initiation sites are determined by auxin 
accumulation in "convergence" points, which lie in a rim 
of epidermal cells around the ground meristem tissue [191 
27, 28, 29 . These single cells with high auxin concentration 
polarize towards the ground meristem and locally trans- 
port their auxin into a cell in the ground meristem tissue. 
This localized inflow triggers the successive polarization of 
PIN distributions along a strand of cells starting from the 
cell with auxin inflow [19]. The strand of polarized cells 
finally extends up to a previously existing strand of po- 
larized cells, building the pre-pattern for the vascular net- 
work. Starting from the petiole of the leaf primordium the 
polarized cells differentiate then into vascular cells [3U] . In 
particular second order veins in Arabidopsis thaliana ex- 
hibit PIN polarization in opposite directions starting from 
a single bipolar cell, which lies in the ground meristem 
below the auxin convergence point in the epidermal layer 
p"9] . This yet unresolved behavior gives rise to the forma- 
tion of closed vein loops when both oppositely polarized 
strands connect to already formed veins. 

A resolution on the origin of bipolar cells is postu- 
lated by examination of a minimal canalization model for 
the polarization of PIN distribution due to auxin supply 
in a one-dimensional strand of cells. Performing a non- 
linear analysis of the model reveals for each single cell 
two uniform stable states considering polarization. One 
resting state, where efflux facilitators are symmetrically 
distributed within the cell, and one polar state character- 
ized by a constant net transport of auxin due to a polar 
localization of PIN proteins. The model predicts auxin 
triggered polarization pulses and fronts as a consequence 
of a dynamic rearrangement of PIN efflux facilitators to- 
wards the polar state. Cells with continuous auxin supply 
can be in a dynamic bipolar state, from which polariza- 
tion fronts travel to both ends of a strand of cells. The 



role of the underlying kinetic processes becomes explicit 
in the course of the non-linear analysis of the polarization. 
An analytic solution result in quantitative predictions on 
the pulse's and front's auxin amplitude depending on the 
kinetic parameters, establishing a basis for detailed exper- 
imental determination. 



2 Model 

To describe how auxin polarizes the distribution of PIN 
during vein formation we focus on a one-dimensional strand 
of cells, see fig. [11(B), assuming that there is no net auxin 
flow perpendicular to the direction of polarization. The 
strand is subdivided into cells of length i numbered by 
n. Every cell is characterized by a single auxin concen- 
tration A(n) and the concentration of membrane bound 
PIN proteins. We distinguish between PIN proteins in the 
bulk of the cell Pb(n) and adsorbed to the cell membrane 
either on the right hand side or on the left hand side of 
the cell, P r (n), and Pi(n), respectively. The auxin con- 
centration per cell changes due to synthesis and degra- 
dation with rates sa and d,A- Furthermore, the amount 
of auxin changes due to a net flow to neighboring cells 
facilitated by PIN proteins embedded in the cell mem- 
branes of the corresponding cell-cell interface. The net 
flow from cell n to cell n + 1 is, hence, given by J(n) = 
eA[A(n)P r (n) — A(n+ l)Pi(n + 1)], where ca denotes the 
transport efficiency rate across the cell-cell interface. The 
full auxin dynamics is then described by, 



dt 



A(n) = s A - d A A(n) - - \J(n) - J(n - 1)] . (1) 



Stating these dynamics for auxin we assume that auxin 
transport is dominantly from one inner cell to the other. 
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Auxin is known to accumulate in cell interiors 0[3T| since 
auxin is a lipophilic weak acid which easily enters cells 
as undissociated acid, its prevailing form at the pH of 5 
in extracellular space. In the interior of plant cells ion 
pumps keep the pH at 7 leading to the ionization of auxin 
and building a concentration gradient that accumulates 
auxin inside cells. As lipid membranes are impenetrable 
for the now charged molecule, auxin is trapped in the 
cell's interior rendering efflux facilitator necessary. Being 
exported by PIN proteins auxin is free to diffuse in the 
extracellular space, the apoplast. However, the distance 
between neighboring cells is so small that almost all auxin 
molecules have entered any cell again within one millisec- 
ond^] Furthermore, auxin can be assumed to be approxi- 
mately uniformly distributed within typical cells because 
the diffusion of the small molecule auxin is very large. 
In aqueous solutions D — 670/im 2 /s has been measured 
|34j which has been confirmed by indirect measurements 
of the diffusion constant in auxin transport experiments 
[35 . Distinct auxin gradient would therefore only arise 
for large plant cells of about lOO^m or larger. It is not en- 
tirely clear that this reasoning holds under the condition 
of highly effective auxin transport |36j , however, the time 
scale of auxin transport by mere diffusion through a cell 
of typical 50/im is 4s, faster than other process contribut- 
ing to the polarization of PIN, substantiating the neglect 
of auxin gradients within cells during PIN polarization. 
Hence, we can approximate auxin flow to be dominated by 
cell-to-cell transport. Synthesis and degradation of auxin 
are considered since they take place on fast time scales 
as the majority of auxin is stored in its conjugated form 
inside the cell which is readily hydrolyzed in less than 
seconds [37] ■ This is in contrast to the production and 
degradation of PIN proteins, which takes place on much 
larger time scales of several minutes. We therefore model 
the total amount of PIN proteins Ptot to be constant per 
cell, yielding the following equality for the number of free 
PIN proteins in the bulk P b (n) = P tot - P r (n) - Pi(n). 
Hence, we only consider the dynamics of PIN proteins em- 
bedded in a cell membrane. Their concentration changes 
first of all by a basal adsorption rate sp and a basal des- 
orption rate dp. Additionally, the net auxin flow over a 
cell-cell interface is modeled to feed back onto the amount 
of PIN proteins favoring the flow direction. This is cast 
in an enhanced attachment or equally a decreased des- 
orption rate: gpJ 2 {n)9(J(n)) as proposed by canalization 
models [TMIMI7JI21Elll3aiffl2B] . By imposing the Heav- 
iside step function 9, the feedback reacts to positive net 



flow only. For the following analysis the feedback is pro- 
portional to the square of the net auxin flow as stated 
above, in the discussion we explain that any exponent 
larger than one yields analogous results. Different feed- 
back mechanisms proposed recently P3, 41,42,321 are also 
compared to our approach in the discussion. Incorporat- 
ing the positive feedback on auxin flow, the PIN dynamics 
are given by, 
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Except for the non-linear feedback term, we assumed through- 
out the model setup linear relationships as a first order ex- 
pansion to the yet elusive detailed underlying dynamics. 
We refer to the discussion for an analysis of model mod- 
ifications confirming the robustness of our assumptions. 
In contrast to many existing canalization models [16,17, 
38,39,40,26 we account for the detailed PIN cycling by 
endosomes similar to Ref. |24j . however, we discard the 
explicit dynamics of a putative auxin synthesizer used in 
that work. 

Up to now our model involves six kinetic rates, how- 
ever, rescaling the concentration of auxin a = A/A eq , 
where A cq — SA/dA, and PIN proteins, p r j = P r ,i/Ptot, as 
well as time r = tdp reveals that only four independent, 
dimensionless parameters govern the behavior of auxin 
and PIN dynamics, namely S a — dA/dp, a p = sp/dp, 
1 P = gpAlqPtoJdp, and e a = e A Pt ot /£dp. Quantitative 
knowledge of the kinetic rates is very sparse. Half-life 
measurements of auxin yield estimates for its degrada- 
tion rate, d A = 2 ■ 1(T 4 - 2 • lCT 5 l/s [HJ, which is how- 
ever strongly affected by environmental conditions such as 
light, wind, and temperature. Permeability measurements 
45 , 32J of PIN assisted auxin anion transport are found 
to be eA-Ptot = lA/im/s. For the other kinetic rates no 
experimental estimates are available to best of our knowl- 
edge although various rates have been assumed in simu- 
lations. This limited knowledge of the kinetic rates un- 
derlying auxin and PIN dynamics demonstrates how de- 
sired an intuition of their relation and role is, which can 
be obtained from mathematical analysis, opening up new 
approaches for experimental measurements. 



1 The importance of extracellular diffusion of auxin can be 
assessed by estimating the residence time of auxin in extracel- 
lular space. Assuming an auxin molecule diffuses with diffusion 
constant D — Q7fim 2 / s 32 in a typical cell-cell interface of 
0.5/im. If it comes close to either of the cells the molecule may 
reenter. Taking into account that not all extracellular auxin 
molecules are protonated and hence able to penetrate the mem- 
brane we assume the probability to enter a cell to be of 10% 
[31IIB] . Considering these assumptions already 97% of all auxin 
molecules have reentered any cell including the one they were 
delivered from after a time duration of one millisecond 133 . 



3 Results 

3.1 Observations from numerics 

During vein formation auxin supplied from the outer epi- 
dermis enters a single cell initiating the polarization of a 
cell strand [Tj5] . We simulated this scenario by integrating 
the microscopic equations (l|3 



kinds of auxin supply, see fig 



) numerically for different 
2] Starting from a strand 
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Fig. 2. Spatial trajectory of polarization pulse and fronts. 
Auxin (blue) and PIN concentration on the left (green) and 
right (orange) hand side of a cell are displayed along a strand 
of cells numbered by their distance in cells from the site of 
initiation. (A) A short, initial supply of auxin to a single cell 
yields a single polarization pulse, while (B) a continuous, high 
inflow of auxin from one end (left) yields a polarization front. 
(C) Continuous supply at a center cell results in two opposite 
polarization fronts originating from a single bipolar cell. Tra- 
jectories arise from the numerical integration of our mathemat- 
ical model for parameters values d^/dp = 0.2, sp /dp — 0.2, 



eAPtot/ldp = 2, and g P At 
auxin supply of A = 20 A; 



P 2 

!*■ tot 



/d T 



8. Initial and continuous 



of cell with evenly distributed PIN proteins and an equi- 
librated amount of auxin, we shortly applied auxin by in- 
creasing the initial auxin concentration in a single cell. 
This triggers a polarization pulse which shortly polarizes 
the PIN proteins in each cell before all cells relax back to 
their initial non-polar state, fig. 2] (A) . Subsequent pulses 
can only be excited when the cells are almost relaxed back 



to their non-polar state therefore a certain lag time is re- 
quired (data not shown) . If the auxin is supplied contin- 
uously by keeping the amount of auxin in a single cell at 
high level a polarization front forms. The front causes all 
cells which it passed to become permanently polarized, 
fig. [2] (B). If a cell in the center of a strand of cells is con- 
tinuously supplied with auxin two fronts arise traveling to 
opposite directions along the strand, fig.[2](C). The latter 
two observations resemble those from vein formation [19] . 
As the polarization pulse and front bear a lot of char- 
acteristics in common, deriving analytic solutions for the 
first gives also quantitative insight into the second. In the 
following our non-linear analysis explains the formation of 
a polarization pulse and front. Identifying the role of the 
underlying kinetic processes by exemplarily solving the 
polarization pulse we derive quantitative results for pulse 
and front properties. Observations in fig. [2] (A) and (B) 
indicate that changes in concentration of PIN proteins on 
the left, Pi, facing adverse to the direction of transport, are 
very small. We therefore assumed in our following analysis 
dPi(n)/dt = 0, i.e., considering the stationary state value 
Pi(n) = <7 P (1 — P r (n))/(l + (jp) for a polarization traveling 
to the right. 



3.2 Static state of a single cell 

Assuming a uniform state for a whole strand of cells each 
single cell itself has two stable and one unstable equilib- 
rium state, as shown in fig. [3] (A) . The first stable fixed 
point at 

RES 



1. 



(4) 



RES 



P 



l + 2a v 



is a resting state, where PIN proteins are evenly distributed 
and no net auxin flow occurs. For parameters beyond 
lp > 4(1 + 2cr p ) two further crossings of the nullclines 
dA/dt — and dP r /dt — occur in a saddle-node bifur- 
cation, a pair of one unstable and one stable fixed point 
at 
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respectively. At the second stable fixed point the PIN dis- 
tribution is polar as P r outnumbers Pi by at least P tot /2, 
yielding a constant net flow of auxin to the right. The rest- 
ing state originates from synthesis and degradation terms 
in eqs. (IK I, while the second pair of fixed points arises 
due to the feedback. As both the resting and the polar 
state are linearly stable a uniform set of cells decays into 
one of them depending on the cells initial state. Left of 
the unstable manifold embedding the unstable fixed point, 
depicted as dashed line in fig. [3] (A), all cells relax to the 
resting state while right to this separatrix all states decay 
to the polar fixed point. This is true for a homogeneous 
set of cells, however, in a spatially inhomogeneous system 
complex scenarios such as waves and fronts arise |46 p 47 ) . 
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Fig. 3. Illustration of the non-linear characteristics of auxin induced polarization. (A) Trajectory of a polarization pulse (grey) 
and front (black) in the nullcline graph obtained from integrating the microscopic equations. Each single cell has two linearly 
stable fixed points (solid green circle), a resting state with symmetric PIN distribution (RES) and a polar state with constant 
auxin flux (POL). Their areas of attraction are separated by an unstable manifold (green dashed line) embedding the unstable 
fixed point (open green circle). Parameters as in fig. [2] (B) Heuristic mapping of the numerical polarization pulse along a strand 
of cells. Auxin is piled up in a cell with even distribution of PIN on both cell membranes at the front of the pulse. Due to the 
slow attachment rates of PIN proteins gp, sp/dp the cell is still non-polar. Increasing the attachment rates due to the positive 
feedback of auxin flow reduces the auxin amplitude. Growing cell-to-cell transport efficiencies eA increase the auxin amplitude 
as more auxin reaches the peak per time step. 
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3.3 Dynamic transition 

Due to spatial inhomogeneities passed on along a strand 
of cells from a cell with auxin supply, the state of a cell 
changes over time as a polarization pulse or front trav- 
els through. The trajectory of states of a single cell in 
time, fig. M (A), maps onto the trajectory of a polariza- 
tion pulse over a strand of cells in space, see fig. [3] (B), 
allowing for a heuristic interpretation of the polarization 
dynamics. To investigate these dynamics we performed 
the continuum limit of the microscopic equations (I1J) , (|2J) , 
see Supplemental I for details. This gives rise to a set of 
partial differential equations, which describe the change 
of a and p r at one point in space over time, as shown in 
eqs. ( 6|7 l above. This time evolution depends on reaction 
terms, including only auxin and PIN concentration at the 
same specific point in space, and gradient terms which ac- 
count for the influence of neighboring sites. The reaction 
terms cause the system to relax to its stable fixed points 
as described in the previous section. The gradient terms, 
however, drive the system along its pulse trajectory. 

The following explains the dynamic transition and the 
role of reaction and gradient terms in it starting from a 
cell in the non-polar resting state (RES). If the gradient in 
auxin and PIN to the neighboring cell is large enough the 
cell is forced out of its stable resting state to larger values 
of auxin entering the domain of attraction of the polar 



state (POL). If a neighboring cell has accumulated more 
auxin and has a higher amount of PIN facing the direction 
of the polarization pulse, p r , it is very effective in trans- 
porting auxin onwards into a cell, raising the auxin content 
well above the equilibrium value. As the now auxin sup- 
plied cell has itself more auxin to transport onwards, the 
net flux increases starting off the positive feedback which 
results in PIN polarization. The then fully polarized cell 
is very efficient in moving its excess auxin onwards, finally 
decreasing its auxin content towards the polar stable state. 
The neighboring cell that has been polar and transport- 
ing on auxin for a bit longer has less auxin, reversing the 
direction of the auxin gradient. If this auxin gradient is 
large enough it drives the cell past its polar state into the 
domain of attraction of the resting state. Hence, if the 
amount of auxin in the neighboring cell is very low, the 
auxin supply breaks down and with it the onward flux of 
auxin. The positive feedback is decreasing and with it po- 
larization towards the non-polar resting state. Subsequent 
polarization pulses can only be triggered if the polariza- 
tion already saturated down to almost resting state values, 
otherwise the remaining polarization would just transport 
the applied auxin onwards before the positive feedback 
can build up an enhanced polarization of PIN proteins. 
The phenomenon that a system has to relax back to its 
resting state before a new pulse can be excited is denoted 
refractory phase in excitable media. 
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We find that the gradient terms in the auxin and PIN 
dynamics have unequal analytic structures that lead to 
their different functions. In the PIN dynamics the squared 
gradient increases spatial inhomogeneities in PIN distri- 
bution by augmenting p r up to saturation. In the auxin 
dynamics the signs of auxin and PIN gradients decide the 
direction of change in auxin concentration. In front of an 
excitation pulse the gradient induces a growth of auxin 
content, while it decreases the amount of auxin at the 
rear of the pulse. The magnitude of the auxin gradient 
also decides between the formation of a polarization pulse 
or front. If the auxin content in the neighboring cell is 
only slightly smaller due to continuous or only slowly vary- 
ing auxin supply a constant flux of auxin through polar 
cells is established, a polarization front forms. If, however, 
the auxin gradient is large and the auxin inflow decreases 
drastically the feedback breaks down and the cells in the 
polarization pulse relax back to their resting state. In sum- 
mary the polarization of PIN distribution by auxin flow is 
initially a bistable system that behaves like an excitable 
medium depending on the amount and continuity of auxin 
supply. 



3.4 Analytical results 



Identifying polarization dynamics as an excitable medium 
enables us to go beyond numerical integration of the mi- 
croscopic equations and analytically compute the auxin 
amplitude of a polarization pulse or front. To this end we 
employed singular perturbation theory 49, 50J on a polar- 
ization pulse. Here we explain the outline of the calcula- 
tion, a detailed derivation is provided in Supplemental I. 
The whole polarization pulse can be subdivided into four 
regions as shown in fig. [2] (A) and fig. [3] (A), first the 
front and back where the auxin concentration is nearly 
constant and only the number of PIN proteins changes 
significantly. Second the excited and the refractory do- 
main during which the efflux facilitator concentration fol- 
lows approximately the nullcline and hence only changes 
according to the nullclines' variation with the auxin con- 
centration while the auxin concentration itself varies pro- 
foundly. Therefore, all four regions are governed to good 
approximation by just a single non-linear equation, the 
one of the PIN protein dynamics or the one of auxin dy- 
namics, respectively. Unfortunately, even each single con- 
tinuum equation is not analytically solvable. The model 
equations ([6j [7]) are therefore linearized around the stable 
nullclines yielding two sets of equations, one left to the 
unstable manifold and one right to it. By this lineariza- 
tion we overestimated the PIN dynamics close to the un- 
stable manifold which leads to smaller auxin amplitudes 
than those resulting from the integration of the micro- 



(A) 



scopic equations ( 1 
by imposing a trave 



3|). An algebraic solution is obtained 
vt), 



ing wave Ansatz A(x — vt) , P r (x~ vt) 
and solving all four equations under the condition of dif- 
ferential continuity at their intersections. The calculation 
yields a closed expression for the auxin amplitude which 
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Fig. 4. Results for the auxin amplitude A max . The auxin am- 
plitude A max varies during a polarization pulse or front with 
its independent kinetic parameters (A) transport efficiency of 
auxin e a = eAPtot/£dp, (B) basal attachment rate a v = sp /dp 
and, (B) enhanced attachment rate 7 P = gpA^^P^/dp. Shown 
are results from the numerical integration of the microscopic 
equations and the analytic expression multiplied by an overall 
factor of 2.3. We considered <j p — sp/dp <C 1 such that only 
less than a third of all PIN proteins occupy each membrane in 
the resting state. e p = eAPtot/£dp > 1 as auxin permeability 
e^Ptot = lAfim/s [431152"] is roughly larger than endosome cy- 
cling by active transport along a cell's cytoskeleton lsp [48] . 
We suggest 7 P = gpA^. q P^ ot /dp > 1 as protein and auxin num- 
bers might be very large. Finally, assuming literature values 
of dA [IU, endosome cycling, and taking cell length of tens 
of jj,m, we used d^/dp = 0.2 to compare our results. Each 
graph shows the variation of a single parameter, while the re- 
maining are kept constant at d^/dp = 0.2, e^Ptot/^dp — 10, 
g P A 2 e<1 P? ot /dp = 12, and s P /d P = 0.2. 
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captures the role of the underlying kinetic parameters: 
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where we abbreviated p = 4(1 + 2a p )/'j p , j3 = 16- 

and 4> = tan^-S^ 2 , ^2 14 - 3 3 /3 4 /). The analytic re- 
sult is compared to the numeric integration of the micro- 
scopic equations 



(1 



3) over a broad range of the kinetic 



parameters in fig. 4 Due to the linearization of the equa- 
tions the algebraic amplitudes are too low. Fitting the an- 
alytic to the numeric results yields an overall factor 2.3. 
This constant factor does not depend on the specific pa- 
rameter range. Considering the amount of approximations 
that entered the calculation, the analytic result captures 
very well the dependence on the kinetic parameters over 
orders of magnitudes. To obtain insight in how the magni- 
tude of the kinetic parameters determines the amplitude, 
we simplified eq. ([8]) further. Following our considerations 
on the size of the kinetic parameters in fig. |4j /3 is smaller 
than one. Hence, expanding for small /3 gives, 
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dp 



(9) 



As the amplitude arises due to auxin inflow from neigh- 
boring cells A max increases with the rescaled transport 
efficiency eAPtot I tdp. Auxin is accumulated in a cell until 
the cell reaches its almost fully polarized state and can 
efficiently transport auxin onwards, see fig. 3KB). Hence, 
accumulation time and auxin amplitude A^^^ inversely 
depend on the basal PIN cycling rate sp/dp and the 
enhanced attachment rate gpP^ ot Al q /dp. Synthesis and 
degradation of auxin do not contribute to the amplitude 
in this first order approximation as their impact on the 
transport is very small. Note, that the amplitude like all 
other pulse and front characteristics is independent of the 
amount of supplied auxin, a general property of excitable 
media. Applying the above simplification to the resulting 
analytic expressions for the velocity of the auxin pulse, see 
Supplemental eq. (S7), (S8), yields, 



v ex e A P u 



(10) 



This result reveals that the pulse of front velocity is dom- 
inated by membrane permeability of auxin facilitated by 
PIN proteins which is just the product of transport effi- 
ciency and number of PIN proteins per cell. The prefactor 
in eq. (10) u/e^Ptot can be estimated from comparison 
with numeric integration of the microscopic equations as 
shown in fig. [I] to be in the range of 6 • 10 -5 — 6 ■ 10~ 4 
depending on the remaining kinetic parameters. 



4 Discussion 

We have shown that some prominent aspects of auxin and 
PIN dynamics can be inferred from a simple mathematical 



model. Our model predicts the transition of ground meris- 
tem cells from a non-polarized stable state to a polarized 
stable state of constant auxin flow. This development oc- 
curs by a traveling wave front triggered by a continuous 
inflow of auxin from the outer epidermal layer, in accor- 
dance with experimental observations [15] . Each cell is a 
bistable excitable medium. Excitations from one state to 
the other can be induced by supply of auxin and crucially 
depend on the spatial gradients in auxin and polar PIN 
concentrations between cells. The amplitude of auxin in 
the wave front and the polarization of the stable states is 
cast in analytical expressions concordant with numerical 
integration of our microscopic equations. 

The microscopic equations underlying our model for 
auxin and PIN dynamics are defined such that all rele- 
vant biological processes are included while a minimum of 
assumptions on their actual kinetics entered. To this end 
only linear synthesis, degradation and transport etc. is 
considered as a first order approximation of any kind of 
kinetics. However, the model is robust against alteration 
of the linear relationships, as is illustrated in Supplemen- 
tal II. For example, extending the cell to cell auxin trans- 
port to account for Michaelis-Menten kinetics preserves 
the form of the nullclines and the dynamics of the wave. 
In our model only the feedback of auxin flow on the attach- 
ment of PIN proteins enters non-linearly. A linear growth 
of the enhanced attachment rate with the auxin flow can- 
not lead to a propagating front, as such a model does not 
exhibit two stable fixed points. Only auxin flow exponents 
higher than one show these properties. However, the exact 
value of the exponent does not affect the form and dynam- 
ics of the traveling wave, again confirming the robustness 
of our assumptions. 

Recently, models for auxin and PIN dynamics were 
developed proposing that the auxin concentration in the 
neighboring cell feeds back onto an enhanced attachment 
rate of PIN proteins 23,41,42,43] in contrast to canaliza- 
tion models, where the net auxin flux governs the feed- 
back. These concentration driven models exhibit a static 
state of spatially ordered auxin maxima with PIN proteins 
polarized towards these auxin maxima |23j . This behavior 
arises as concentration driven feedback changes the non- 
linear character of auxin and PIN dynamics. These models 
generally exhibit only a single stable, resting fixed point 
51,52J. The polarization due to auxin supply observed in 
these models E 42J arises due to an evolved relaxation into 
the stable, resting state. Hence, the polarization is only 
temporary and, for instance, the amount of polarization 
and the velocity of the polarization front depend crucially 
on the amount of auxin supply. This is in contrast to our 
minimal canalization model, where all polarization char- 
acteristics are only governed by the kinetic parameters. 
The amount of auxin supply in an excitable medium only 
regulates if a pulse or front is excited or not. Hence, these 
qualitative differences may help to distinguish between the 
different models experimentally. 

The role of all kinetic processes during the dynamic 
rearrangement of PIN and auxin in cells becomes explicit 
when examining the very front of the polarization in a mi- 
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croscopic scenario as illustrated in fig. [31(B). The almost 
fully polarized cell at the peak of the front carries a lot of 
auxin molecules that are invading the next yet non-polar 
cell in the direction of polarization with a rate mainly 
governed by the cell-to-cell transport efficiency e^. To suc- 
cessfully transfer the accumulating auxin onwards the PIN 
proteins in the yet non-polar cell have to rearrange to fa- 
cilitate directed transport. However, the endosome cycling 
sp, gp by which the membrane bound PIN proteins reach 
the cell membrane is very slow. Hence, the attachment 
rate of the efflux facilitators forms a bottleneck that piles 
up more and more auxin in a cell, that is slowly increasing 
the amount of PIN proteins facing the direction of trans- 
port. Heuristically, an auxin pulse forms due to a traffic 
jam caused by the slow cycling of the efflux facilitators, as 
given by eq. (pM . As an equilibrium concentration of PIN 
proteins is always embedded in each membrane ready to 
transport auxin, the magnitude of the velocity of a po- 
larization front or pulse is set by the cell-to-cell transport 
efficiency, see eq. (10). The other kinetic parameters only 
slightly modulate the velocity. The PIN attachment rates, 
gp, sp, and detachment rate, dp, on the other hand de- 
termine the number of PIN proteins accumulated at the 
membrane in the polar stable state, see eq. ^j. The po- 
larization grows with the enhanced attachment rate, the 
strength of the feedback, gp. On the contrary, basal PIN 
cycling sp/dp intensifies the competition between oppos- 
ing membranes decreasing the amount of PIN proteins in 
the direction of polarization. 

The result of our analytic expressions for the PIN con- 
centration in the polar state with constant auxin flux in 
eq. (pi) and the auxin amplitude at the very head of the po- 
larization front eq. ([8]) enable estimates of the underlying 
kinetic rates by identifying and measuring these observ- 
ables in future experiments. Existing experimental results 
by Scarpella et al. [19) permit an estimate of the polar- 
ization front velocity in the range v = 10 -4 — 10 _3 /zm/s, 
in accordance with our estimate for the velocity v — 8 • 
1CP 5 — 8 • 10~ 4 /im/s, resulting from the fitted pre-factor 
in eq. \1Q\ and the literature value of auxin permeability 
e a Ptr^t |4*5ll3"2"] . An quantitative estimate of PIN polarity 
from the same existing data is to best of our knowledge yet 
unfeasible as a reference for the protein number is absent. 
This could be overcome by new experiments, which could 
also aim at the auxin kinetic rates. Unlike PIN which is 
readily GFP tagged, auxin is not directly detectable and 
quantification of its amount can only occur via indirect 
methods. Recently, measurements of deuterated auxin im- 
proved [5 3) making experiments with exogenously applied 
auxin conceivable. In such setups one should, however, 
keep in mind that exogenous auxin mixes with endoge- 
nous, non-labelled auxin, decreasing the observed ampli- 
tude. The position of the auxin peak can easily be located 
as it should be accumulated in front of those cells with 
the largest amount of PIN proteins at the corresponding 
membrane. Measurements of the amount of PIN proteins 
in polarized cells via GFP tagging could not only disclose 
the basal endosome cycling rate but also the magnitude 
of the feedback between auxin flow and PIN dynamics. 



The occurrence of bipolar cells has stimulated pre- 
vious theoretical models introducing a hypothetical new 
molecule [54 or moving auxin sources [38] , Our model, 
however, readily predicts the occurrence of bipolar cells 
along a one-dimensional strand of cells at the site of con- 
tinuous auxin inflow. These cells show a high concen- 
tration of PIN proteins on either membrane, a balanced 
outcome of the competition for PIN between both mem- 
branes. This state is not a statically stable but dynami- 
cally driven by the supply of auxin. Transferring this ob- 
servation to the two-dimensional layer of ground meristem 
cells during vein initiation may explain the bipolar cells 
observed experimentally [19,27,28,23]. In two dimensions 
several membranes can compete, yielding also triple polar 
cells or theoretically higher orders of polarity. However, 
in biological cells not all cell membranes may have the 
same number of PIN proteins to enter the competition. 
Those with fewer initial PIN proteins will become the site 
of PIN drain, decreasing the number of successful polar- 
ization fronts. Experiments indicate that mechanical cues 
might favor certain membranes |55| paving the way of a 
polarization front and hence the position of veins and vein 
loops. 

In summary we analyzed a canalization model to ex- 
plain how auxin and PIN dynamics polarize the distribu- 
tion of PIN proteins during early vein initiation in the 
ground meristem. Each cell is found to be bistable consid- 
ering their PIN distribution in the membrane. Polarization 
occurs in a traveling front as auxin is supplied to a single 
cell mimicking the auxin inflow from the outer epider- 
mal cell layer. The driver in this transition is the spatial 
gradient of auxin and polar PIN concentration between 
cells. The key idea in this polarization is the positive feed- 
back between hormone auxin and its own efflux facilitator 
PIN. An idea that might be inspiring in other developmen- 
tal processes in animals where tissue is polarized as, for 
example, in planar cell polarity |56j . As the polarization 
front can travel in opposite directions from its initiation 
site, a bipolar cell, the up to now puzzling occurrence of 
closed vein loops can be resolved. Furthermore, our non- 
linear analysis enables the calculation of exact analytical 
expression for the polarization front. Therefore, our new 
quantitative predictions for the PIN polarization driven 
by auxin flow establishes a basis to determine the kinetic 
parameters underlying the transport of auxin and may 
therewith have far-reaching impacts on the understand- 
ing of the developmental processes and their differences in 
plant species, to perceive the fundamental patterns of leaf 
veins or phyllotaxis and learn how environmental condi- 
tions alter these. 



Acknowledgements 

The authors thank Eric M. Kramer for fruitful discussions. 
Financial support of the German Excellence Initiative via 
the program "Nanosystems Initiative Munich (NIM)" and 
of the LMUinnovativ project "Functional Nanosystems 
(FuNS)" is gratefully acknowledged. K.A. acknowledges 
funding by the Studienstiftung des deutschen Volkes. 



K. Alim and E. Frey: Auxin-induced polar distribution of PIN proteins 



References 

1. J. Petrasek, J. Friml, Development 136(16), 2675 (2009) 

2. J. Friml, J. Wisniewska, E. Benkova, K. Mendgen, 
K. Palme, Nature 415(6873), 806 (2002) 

3. A. Marchant, J. Kargul, S.T. May, P. Muller, A. Delbarre, 

C. Perrot-Rechenmann, M.J. Bennett, EMBO J 18(8), 
2066 (1999) 

4. D. Reinhardt, E.R. Pesce, P. Stieger, T. Mandel, K. Bal- 
tensperger, M. Bennett, J. Traas, J. Friml, C. Kuhlemeier, 
Nature 426(6964), 255 (2003) 

5. L.E. Sieburth, Plant Physiol. 121(4), 1179 (1999) 

6. F.A.F.C. Went, Naturwissenschaften 21, 1 (1933) 

7. P.H. Rubery, A.R. Sheldrake, Nature - New Biol. 
244(139), 285 (1973) 

8. M.H.M. Goldsmith, T.H. Goldsmith, M.H. Martin, Proc. 
Natl. Acad. Sci. USA 78(2), 976 (1981) 

9. G.J. Mitchison, Proc. Roy. Soc. Lond. B Bio. 209(1177), 
489 (1980) 

10. L. Galweiler, C.H. Guan, A. Muller, E. Wisman, K. Mend- 
gen, A. Yephremov, K. Palme, Science 282(5397), 2226 
(1998) 

11. K. Okada, J. Ueda, M.K. Komaki, C.J. Bell, Y. Shimura, 
Plant Cell 3(7), 677 (1991) 

12. E. Benkova, M. Michniewicz, M. Sauer, T. Teichmann, 

D. Seifertova, G. Jiirgens, J. Friml, Cell 115(5), 591 (2003) 

13. J. Friml, A. Vieten, M. Sauer, D. Weijers, H. Schwarz, 
T. Hamann, R. Offringa, G. Jiirgens, Nature 426(6963), 
147 (2003) 

14. J. Petrasek, J. Mravec, R. Bouchard, J.J. Blakeslee, 
M. Abas, D. Seifertova, J. Wisniewska, Z. Tadele, 
M. Kubes, M. Covanova et al., Science 312(5775), 914 
(2006) 

15. T. Sachs, Ann. Bot. - London 33(130), 263 (1969) 

16. G.J. Mitchison, Proc. Roy. Soc. Lond. B Bio. 207(1166), 
79 (1980) 

17. G.J. Mitchison, Philos. Trans. Roy. Soc. B 295(1078), 461 
(1981) 

18. M. Sauer, J. Balla, C. Luschnig, J. Wisniewska, V. Reinohl, 
J. Friml, E. Benkova, Genes Dev 20(20), 2902 (2006) 

19. E. Scarpella, D. Marcos, J. Friml, T. Berleth, Genes Dev 
20(8), 1015 (2006) 

20. T. Paciorek, E. Zazimalova, N. Ruthardt, J. Petrasek, 
Y.D. Stierhof, J. Kleine-Vehn, D.A. Morris, N. Emans, 
G. Jiirgens, N. Geldner et al., Nature 435(7046), 1251 
(2005) 

21. P. Dhonukshe, I. Grigoriev, R. Fischer, M. Tominaga, D.G. 
Robinson, J. Hasek, T. Paciorek, J. Petrasek, D. Seifer- 
tova, R. Tejos et al., P Natl Acad Sci USA 105(11), 4489 
(2008) 

22. R. Benjamins, B. Scheres, Annual Rev Plant Biol 59, 443 
(2008) 

23. H. Jonsson, M.G. Heisler, B.E. Shapiro, E.M. Meyerowitz, 

E. Mjolsness, Proc. Natl. Acad. Sci. USA 103(5), 1633 
(2006) 

24. F.G. Feugier, A. Mochizuki, Y. Iwasa, J. Theor. Biol. 
236(4), 366 (2005) 

25. T. Berleth, E. Scarpella, P. Prusinkiewicz, Trends Plant 
Sci 12(4), 151 (2007) 

26. R.S. Smith, E.M. Bayer, Plant Cell Environ 32(9), 1258 
(2009) 

27. C.L. Wenzel, M. Schuetz, Q. Yu, J. Mattsson, Plant J 
49(3), 387 (2007) 



28. E.M. Bayer, R.S. Smith, T. Mandel, N. Nakayama, 
M. Sauer, P. Prusinkiewicz, C. Kuhlemeier, Genes Dev 
23(3), 373 (2009) 

29. D. Koenig, E. Bayer, J. Kang, C. Kuhlemeier, N. Sinha, 
Development 136(17), 2997 (2009) 

30. M.G. Sawchuk, P. Head, T.J. Donner, E. Scarpella, New 
Phytol 176(3), 560 (2007) 

31. J.A. Raven, New Phytol. 74(2), 163 (1975) 

32. R. Swarup, E.M. Kramer, P. Perry, K. Knox, H.M.O. 
Leyser, J. Haseloff, G.T.S. Beemster, R. Bhalerao, M.J. 
Bennett, Nat. Cell. Biol. 7(11), 1057 (2005) 

33. A.S. Usenko, A.G. Zagorodny, Mol. Phys. 61(5), 1213 
(1987) 

34. P. Larsen, Growth substances in higher plants, Vol. 3 
of Modern methods of plant physiology (Springer, Berlin, 
1955) 

35. E. Wangermann, G.J. Mitchison, Plant Cell Environ. 4(2), 
141 (1981) 

36. E.M. Kramer, J Exp. Bot. 59(1), 45 (2007) 

37. S. Pollmann, A. Muller, M. Piotrowski, E.W. Weiler, 
Planta 216(1), 155 (2002) 

38. A.G. Rolland-Lagan, P. Prusinkiewicz, Plant J 44(5), 854 
(2005) 

39. S. Stoma, M. Lucas, J. Chopard, M. Schaedel, J. Traas, 
C. Godin, PLoS Comput Biol 4(10), el000207 (2008) 

40. E.M. Kramer, Trends Plant Sci 14(5), 242 (2009) 

41. R.S. Smith, S. Guyomarc'h, T. Mandel, D. Reinhardt, 
C. Kuhlemeier, P. Prusinkiewicz, Proc. Natl. Acad. Sci. 
USA 103(5), 1301 (2006) 

42. R.M.H. Merks, Y.V. de Peer, D. Inze, G.T.S. Beemster, 
Trends Plant Sci 12(9), 384 (2007) 

43. M. Ibanes, N. Fabregas, J. Chory, A. Cano-Delgado, P. 
Natl. Acad. Sci. USA 106(32), 13630 (2009) 

44. F. Rapparini, YY. Tarn, J.D. Cohen, J. P. Slovin, Plant 
Physiol. 128(4), 1410 (2002) 

45. A. Delbarre, P. Muller, V. Imhoff, J. Guern, Planta 198(4), 
532 (1996) 

46. V.S. Zykov, A.T. Winfree, Simulation of Wave Processes 
in Excitable Media (Manchester Univ. Press, Manchester, 
1987) 

47. A.S. Mikhailov, AY. Loskutov, Foundations of Synergetics 
I: Distributed Active Systems (Springer Verlag, New York, 
1991) 

48. C.L. Howe, W.C. Mobley, J Neurobiol 58(2), 207 (2004) 

49. P. Ortoleva, J. Ross, J. Chem. Phys. 63(8), 3398 (1975) 

50. R.G. Casten, H. Cohen, P.A. Lagerstrom, Q. Appl. Math. 
32(4), 365 (1975) 

51. A. Newell, P. Shipman, Z. Sun, J Theor Biol 251(3), 421 
(2008) 

52. P. Sahlin, B. Soderberg, H. Jonsson, J Theor Biol 258(1), 
60 (2009) 

53. K. Ljung, A.K. Hull, J. Celenza, M. Yamada, M. Estelle, 
J. Normanly, G. Sandberg, Plant Cell 17(4), 1090 (2005) 

54. F. Feugier, Y. Iwasa, J Theor Biol 243(2), 235 (2006) 

55. O. Hamant, M.G. Heisler, H. Jonsson, P. Krupin- 
ski, M. Uyttewaal, P. Bokov, F. Corson, P. Sahlin, 
A. Boudaoud, E.M. Meyerowitz et al., Science 322(5908), 
1650 (2008) 

56. P.A. Lawrence, G. Struhl, J. Casal, Nat Rev Genet 8(7), 
555 (2007) 



10 K. Alim and E. Frey: Auxin-induced polar distribution of PIN proteins 



Supplementary Information: Quantitative Predictions on Auxin-Induced 
Polar Distribution of PIN Proteins during Vein Formation in Leaves 

Karen Alim and Erwin Frey 
Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, 
Ludwig-Maximilians-Universitdt, Theresienstr. 37, D-80333 MUnchen, Germany 

I. AUXIN AMPLITUDE COMPUTATION 



Starting from the microscopic definition of a model for the polarization of PIN distributions 
by auxin flow, see Eq. (1), (2) and (3), an exact analytical expression for the auxin amplitude and 
the velocity of the polarization pulse is derived. To this end a continuum limit is performed and 
singular perturbation theory is employed as described in the following. 

Initially, polarization is defined by three equations governing the dynamics of auxin concen- 
tration A and the amount of PIN efflux facilitators on the membrane on the right P r or on the 
left Pi hand side of the cell. Assuming the symmetry of the system is broken such that polariza- 
tion evolves to the right, the concentration of efflux facilitators on the left hand side membrane of 
every cell does not change significantly with time if s P /(dp + 2s P ) <C 1, we therefore assume 
dPi(n)/dt = 0. As flow proceeds to the right J(n) > Vn, the amount of PIN proteins on the left 
hand side membrane amounts in its stationary state to Pi(n) = Sp(l — P r (n))/(dp + sp). Substi- 
tuting this result in the remaining dynamic equations the system is described by two components 
only, A(n) and P r (n). We derive continuum equations by setting n — > x, n + 1 — > x + £, and 
n — 1 — > x — £. When the wavelength as the length scale on which the pulse evolves is consider- 
ably larger than the cell length £, a time scale separation occurs which makes higher order terms 
negligible small. As observed in the simulations shown in Fig. 2 the wavelength of a single pulse 
is of the order of tens of cells justifying a Taylor expansion in x. To describe the characteristics 
observed in the microscopic equations with continuous equations only zeroth order terms and a 
single second order term in the PIN dynamics are required, 
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In these expressions the degree of nonlinearity is still too high to obtain analytical results by 
use of singular perturbation theory. Therefore, the continuum equations are further simpli- 
fied by approximating them right and left of the nullcline embedding the unstable fixed point, 
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p r( x ) = 2 a (C+2s P ) ~ 2(T P l S 2 P sp) \/ 1 ~ gp££[J) > which to S ood approximation resembles the sep- 
aratrix between the areas of attraction of the stable fixed points. 

Left of the nullcline the reaction terms of the continuum equations, those terms without spatial 
or temporal derivatives, are expanded around the stable nullcline embedding the stable resting state 
(A/A eq = 1, P r /P tot = sp/(dp + 2sp)). The spatial derivative terms are simplified by discarding 
the spatial derivative of auxin whose factor P r (x)/P tot — sp/(dp + 2s p) turns the whole term 
negligible small close to the stable nullcline, resulting in, 
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Right of the nullcline embedding the unstable fixed point the reaction term of the efflux facil- 
itator dynamics is expanded around the polar fixed point (A/A eq = l,P r /P tot — nf?^, 3 *?^ -I 
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1 P X pi )■ Here, the spatial derivative terms contributing to the PIN protein dy- 



namics become negligible small and are therefore discarded, yielding, 
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Based on these sets of equations the amplitude and the velocity of a polarization pulse are calcu- 
lated employing singular perturbation theory [1, 2]. The singular perturbation approach becomes 
applicable to a two component system performing a pulse if the pulse can be separated in different 
regions which comply either of the following restrictions. Either to good approximation one com- 
ponent is constant while the other changes rapidly. Or if both components change simultaneously 
one of the components should follow a nullcline. Then the two coupled nonlinear equations decou- 
ple in each region and only one differential equation remains to be solved. Assuming a traveling 
wave ansatz z = x — vt, where v defines the velocity of the wave, the partial differential equations 
simplify to ordinary differential equations. These differential equations remain to be solved under 
the condition of continuity and differential continuity at the nullcline embedding the unstable fixed 
point which separates the two cases Eqs. (S3, S4) and Eqs. (S5, S6). 

The trajectory of an auxin pulse can be subdivided into four regions, first a wave front and 
back, where the auxin concentration is approximately constant ^4 max ,min, while the PIN protein 
concentration changes rapidly and second an excited and a refractory region during which the PIN 
concentration follows the stable nullclines. Several boundary conditions arise from the requirement 
of continuity and differential continuity. Considering the efflux facilitator dynamics continuity 
requires that PIN concentrations during pulse front and back governed by Eqs. (S4) and (S6) merge 
into the nullclines defining refractory and excited domain P r fi OXA {z — > oo) = -P r ,refrac, P r ,iront{z — > 

-OO) = P r ,excite and P r>b ack(-Z ~> ~00) = Pr.refrac, P r ,back(z -»• 00) = P r , eX cite- As the Separating 

nullcline is crossed during wave front and back additionally continuity and differential continuity 
is compulsory for efflux facilitators at the position of the nullcline embedding the unstable fixed 
point z sep . The auxin concentration evolving during excited Eq. (S5) and refractory region Eq. (S3) 
has to reach the constant auxin concentration of pulse front and back A max , A m i n at distinct points 
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in evolution z = Zf ront and Zback, resulting in the boundary conditions A- e frac( Zfront) = Ai 

^excite v ^front/) ^refrac x ^back/ ^min ^excite v ^back/> O^refrac l^front// 0* O^-excite v *2front//(Z^> anu 

^Aefrac(^back)/^^ = ^excite (zback)/dz. Two of those boundary conditions yield equations that 
solve for the auxin amplitude A max and the pulse velocity v. First the condition of differential 
continuity of the PIN concentration at the separatrix dP r (z sep )/dz results in, 
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where we abbreviated p = 4(d P + 2s P )/ g P Pf oV Second differential continuity of auxin left and 
right of a wave front or back A re frac( Zfront) = -4max = ^excite (-^front) yields, 
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These two equations together result in a quartic equation for the amplitude of the auxin pulse, 
evaluated to the expression given in eq. (8). Reentering this result into one of the defining equations 
above yields an analytic result for the velocity of an auxin pulse. 

II. COMPARISON TO ALTERNATIVE MICROSCOPIC TRANSPORT MODELS 

When defining our microscopic equations Eqs. (1 - 3) we included all contributing processes 
but considered a minimum of assumptions on the kinetics. In this spirit all processes are modeled 
by linear relations as the first order term of any kind of underlying kinetics. The only exception 
is the enhanced attachment of PIN proteins, which is the point where non-linearity enters the 
microscopic model triggering the non-linear effect of a traveling pulse or front. Including further 
non-linearities renders the microscopic models intractable for analytical calculations. However, the 
analysis of our minimal model has revealed the key characteristics for polarization to be two stable 
fixed points accompanied by an unstable manifold. Excitations beyond this unstable manifold 
then lead to the development of a wave pulse or the relaxation to the polar stable fixed point. 
With this knowledge we can assess more evolved microscopic models by comparing their non- 
linear characteristics such as the nullclines to the minimal model. As is shown in the following, 
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FIG. SI: Comparing the polarization pulse for different extended models to our minimal model yields no 
difference in their characteristics confirming the robustness of our linear approximations. Parameters values 
as follows. Minimal model d A /dp = 0.2, sp/dp = 0.2, e A P tot /£dp = 10, and gpA^P^/dp = 12. 
Michaelis-Menten kinetics for active auxin transport d A /dp = 0.2, sp/dp = 0.2, e A P tot /£dpA eq = 8, 
gpP^i/dp = 70, and k A /A eq = 1.6. Feedback to the power of three d A /dp = 0.2, sp/dp = 0.2, 
e A Ptot/ed P A eq = 10, and gpP^A^/dp = 18. 



exemplarily changing linear terms of the minimal model into non-linear terms results only in slight 
changes of the characteristics confirming the robustness of a first order assumption in the minimal 
model. If not explicitly stated otherwise we assumed dPi(n)/dt = to calculate the fixed points. 

A. Michaelis-Menten kinetics for active auxin transport 

One may assume that a Michaelis-Menten mechanism describes the active transport of a sub- 
strate auxin by an enzyme presented by the PIN proteins [3, 4]. Then, in the definition of the 
net auxin flow the number of transported auxin molecules is represented by a Hill function with 
Michaelis-Menten constant k A . For reasons of completeness we state the full set of microscopic 
equations: 
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s A - d A A(n) - e -j [J(n) - J(n - 1)] , 

-d P P r (n) + s P P b (n) + g P J 2 (n)6(J(n))P b (n), 

-dpP t {n) + s P P b {n) + g P J 2 (n - 1)6(-J{n - l))P 6 (n), 
A(n) A(n + 1) 



-PJn 



A(n) + k A A(n + l) + k A 



P/fn + 1 



(S9) 
(S10) 
(Sll) 
(S12) 



These equations yield in accordance with the minimal model three fixed points, one resting state 
at (A/A eq = 1, Pr/Ptot = sp/(dp + 2sp)) and a pair of stable and unstable constant current fixed 



points at (A/A eq = l,P r /P tot = ^±^ T ^^ Jl - (1 + k A ?^g^). The model 
displays very similar nullclines and the same dynamics as the minimal model as exemplified in 
Fig. SI. If the new parameter k A lies outside its range k A < — 1 — \ A u P xj^ n and fc^ > — 1 + 
, '"' s only a single stable fixed point occurs and no polarization can be observed. 



4(dp+2s P ) 



B. Feedback - power of the current 

In the minimal model we take the feedback of auxin flow on the enhanced attachment of PIN 
proteins to enter with a power of two. In general one could assume any kind of power, 

j t A(n) = s A - d A A(n) - e -j [J(n) - J(n - 1)] , (S13) 

j t Pr{n) = -d P P r {n) + s P P b {n) + g P J k {n)6{J{n))P b {n), (S14) 

j t P(n) = -d P P{n) + s P P b (n) + g P J k (n - 1)9(-J(n - l))P b (n), (S15) 

J(n) = A{n)P r {n)-A{n + l)Pi(n + l). (S16) 



Considering only integer powers for simplicity, we find the following. For k = 1 the above 
equations display the resting fixed point at (A/A eq = 1, P r /P tot = Pi/P lo t = s P /{d P + 2sp)) and 
a polar fixed point at (A/A eq = 1, Pi/P tot = sp/g P A eq P tot , P r /P lo t = 1 - (s P + dp)/g P A eq P tot ). 
The resting state is the only fixed point and stable for gpA eq P tot < d P + 2s P , otherwise the polar 
fixed point is stable and the resting state turns unstable. No excited polarization can occur in either 
case. However, if k > 1 the equations display a set of stable and unstable fixed points in addition 
to the resting stable fixed point. For any k > 1 one recovers the dynamics observed for k = 2 in 
the minimal model as exemplified in Fig. SI for the cases k = 3. This observation is in accordance 
with results of Ref. [5], which stated vein patterns for any feedback function obeying to first order 
a higher power than k = 1. For k = 2 we observed that the pair of a stable and an unstable polar 
fixed point occurs only for gpALP^Jdp > 4(1 + 2sp/dp), similar rules apply for k = 3 or higher. 



III. MODEL PARAMETERS 

Our model depends on four dimensionless parameters 5 a = d A /d P , a p = s P /d P , 
1p — fl'P^eq-Ptot/^-P' an d e a = e A P tot /£d P . As quantitative knowledge is very sparse, i.e., 
only d A = 2 • 10~ 4 — 2 • 10~ 5 l/s and e A P to t = 1.4/xm/s are experimentally verified, their 
values have been varied over large ranges within conceptional limits, see table below. As only 
less than a third of all PIN protein can occupy each membrane in the resting state given by 
Pr,l/Ptot — °p/(l + 2<7 P ), <7p is limited to values considerably smaller than one. Furthermore, it 
is reasonable to assume that e A P tot /£d P > 1 as auxin permeability e A P tot = 1-4/um/s [6, 7] is 
roughly larger than endosome cycling by active transport along a cell's cytoskeleton £s P [8]. In 
addition, we suggest g P ALP^ t /d p > 1 as protein and auxin numbers might be very large. Finally, 
assuming literature values of d A [9], endosome cycling, and taking cell length of tens of \xm, we 
took d A /d P < 1. The parameter assumed in our simulation presented in Fig. 5 are summarized in 
the following. 



Parameter 


5 a = d A /d P 


(jp = s P /d p 


lv = dpA^PfJdp 


e a = e A P tot /£d P 


Values 


0.1-1 


0.05 -1 


1-10 


1-100 
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